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SECTION  I 
INTRODUCTION 

By  the  very  nature  of  a  numerical  approach  the  computation 
of  a  flow  field  is  restricted  to  a  finite  region.  Figure  1.  At 
the  outer  boundary  of  this  region  one  must  impose  conditions  which 
insure  that  a  further  continuation  of  the  flow  field  has  the 
properties  desired  at  infinity.  In  a  steady  flow  this  means  that 
the  deviations  from  a  parallel  flow  tend  to  zero  as  one  goes  to 
infinity.  In  an  oscillatory  flow  this  criterion  is  not  sufficient, 
for  it  is  satisfied  by  incoming  as  well  as  by  outgoing  waves.  One 
must  express,  somehow,  that  only  outgoing  waves  are  permissible. 

The  present  discussions  are  carried  out  with  a  view  toward 
the  computation  of  flow  fields  in  the  lower  transonic  region.  The 
flow  in  the  far  field  does  not  show  the  difficulties  which  are 
typical  for  the  transonic  regime;  it  can  be  treated  by  means  of  the 
potential  equation  for  unsteady  flows  linearized  for  the  vicinity 
of  a  parallel  flow  with  the  prescribed  free  stream  Mach  number. 

In  steady  flow,  one  deals  in  essence,  with  .e  Laplace  equation. 

If  oscillations  with  a  fixed  frequency  are  present,  one  ultimately 
has  to  solve  the  Helmholtz  equation. 

The  far  field  conditions  for  steady  transonic  flows  have  been 
studied  by  Klunker  (Ref.  1) .  The  extension  to  oscillatory  flows 
is  due  to  Traci,  Albano,  Farr  and  Cheng  (Ref.  2) .  In  one  regard 
these  formulations  have  an  unexpected  form.  The  conditions  at  the 
far  boundary  express  a  property  of  the  far  field.  It  should, 
therefore,  be  possible  to  formulate  them  in  terms  of  quantities 
pertaining  to  the  far  field,  specifically  in  terms  of  the  potential 
and  its  normal  derivative  at  its  inner  boundary.  These  quantities 
are,  of  course,  identical  with  the  potential  and  the  normal 
derivative  at  the  outer  boundary  of  the  computed  part  of  the  flow 
field.  The  formulations  of  Ref.  1  and  Ref.  2  express  the  far 
field  conditions  in  terms  of  near  field  data.  A  formulation  solely 
in  terms  of  far  field  data  is  indeed  possible.  It  is  not  more 
complicated  than  that  of  Ref.  1  and  2.  We  shall  give  a  rather 
thorough  discussion  of  a  number  of  different  formulations  including 
that  of  Klunker  and  Traci. 


A  clear  understanding  of  the  far  field  conditions  is  particularly 
important  for  oscillatory  flows.  In  a  steady  flow  the  conditions 
at  the  far  boundary  have  only  a  relatively  small  influence.  They 
can  be  taken  into  account  by  iteration.  For  this  purpose  the 
procedures  of  Ref.  1  and  2  are  well  suited.  Actually  they  have 
been  used  only  in  this  sense.  If  one  has  a  flow  with  harmonic 
oscillations,  then  the  far  field  condition  has  a  drastic  effect 
unless  the  frequency  is  very  low.  Take  for  instance  the  field  around 
a  pulsating  sphere  in  air  at  rest  (Helmholtz  equation) .  The  solution 
which  admits  only  outgoing  waves  is  quite  different  from  the  one 
where,  at  a  large  radius,  the  amplitude  of  the  oscillation  is 
required  to  be  zero.  For  such  flows,  it  is  likely  that  one  will 
have  convergence  difficulties  and  under  these  circumstances 
alternative  formulations  of  the  far  field  conditions  may  be  useful. 

The  formulation  of  far  field  conditions  obviously  requires 
that  the  oscillatory  part  of  the  solution  in  the  far  field  can  be 
represented  analytically.  This  is  possible,  only  if  the  governing 
equation  is  rather  simple.  The  Helmholtz  equation,  for  instance, 
has  simple  fundamental  solutions  for  the  far  field.  The  differential 
equations  for  oscillations  superimposed  to  a  flow  with  a  free 
stream  Mach  number  one  and  linearized  for  the  vicinity  of  such  a 
flow  are  much  more  complicated  and  analytic  solutions  are  not 
available.  Such  problems  are  not  included  in  the  present  analysis. 

The  primary  goal  of  this  report  is  the  development  of  the 
basic  concepts.  The  work  is,  of  course,  done  with  a  view  toward 
numerical  implementation;  however,  it  is  still  several  steps 
removed  from  an  actual  program.  An  example  is  given,  but  it  is 
overly  simple.  Its  sole  purpose  is  to  illustrate  the  typical  steps 
which  one  has  to  carry  out.  But  it  also  shows  some  specific 
difficulties  encountered  in  solving  the  Helmholtz  equation.  This  is 
desirable  because  the  experience  gained  with  the  Laplace  equation 
is  not  always  a  safe  intuitive  guide. 

Section  II  starts  with  the  linearized  equation  for  unsteady 
flow,  carries  out  transformations  which  bring  the  problem  into  the 
form  of  the  Helmholtz  equation  and  derives  particular  solutions 
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which  presumably  represent  outgoing  (or  incoming)  waves.  Section 

III  has  a  more  theoretical  nature.  It  shows  that  the  particular 
solutions  derived  in  Section  II  are  a  superposition  of  outgoing 
(or  incoming)  acoustical  waves;  that  is,  it  shows  that  these 
particular  solutions  have  indeed  the  desired  characters.  Section 

IV  uses  Green's  formula  to  discuss  certain  properties  of  the  solutions 
at  the  far  boundary.  The  far  field  conditions  are  fully  formulated 

in  Section  V.  This  includes  a  discussion  of  the  formulations  of 
Ref.  1  and  2.  In  Section  VI  the  principal  steps  needed  in  an 
application  are  discussed  on  the  basis  of  an  example.  It  is  very 
simple  so  that  the  steps  which  usually  are  carried  out  numerically 
can  be  replaced  by  analytical  formulae.  In  these  examples,  one 
recognizes  certain  difficulties  which  may  arise  in  practical  work. 
Sections  VII  and  VIII  contain  a  number  of  observations  concerning 
the  practical  side  of  the  numerical  work. 


SECTION  II 

BASIC  EQUATIONS  AND  PARTICULAR  SOLUTIONS 


The  linearized  equation  for  unsteady  compressible  flow  is 
given  by 

iff  "  f/^tt  (1) 

x,  y,  z  is  a  system  of  Cartesian  coordinates,  t  the  time,  $  the 
velocity  potential,  U  the  freestream  velocity  and  “a ’  the  velocity  of 
sound  (constant  because  of  the  linearization) .  We  shall  write 


M  =  U/a 


One  sets 


where  v  is  a  circular  frequency,  and  x,  y,  z  arise  by  the  Prandtl 
Glauert  coordinate  distortion 

x  '  y  =  /  ,  £  *  0 -tf)  ^ z  (4) 


This  gives 

7  i  j  lcv  Ur 

w  * 


x-  It*-  / 


The  hypothesis 

r'  j  /  tv  U  ~  i 

4  =  $  **/*( dj^u*-*)  (6) 

removes  the  term  <j>  .  The  coordinates  x,  y,  z  are  made  dimensionless 
with  some  characteristic  length  L  (either  the  chord  or  the  half 
chord  of  an  airfoil) . 


y ■  ?A  j 

Then  one  obtains  the  Helmholtz  equation 


4kx  *  Ay  '  At  Y‘¥  -  ° 


with 


olv  L 
a *--  tf 
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Frequently,  one  introduces  in  such  problems  a  reduced  frequency 
to  be  noted  here  by  y 1 

^  -  *L/V 

one  has  .  ^ 

Let  r  and  r'  be  position  vectors  in  the  x,  y,  z  system  and  let 

r^/r-  ?'/ 

A  fundamental  solution  to  be  denoted  by  w(r,  r')  is  obtained  by  the 

coft  r)  *  f(r) 

and  one  obtains 

to(?  r'}  =  - 
or 


(10) 

(11) 


The  expressions  10  and  11  represent  outgoing  and  incoming  waves, 
respectively.  At  least  they  would  represent  outgoing  and  incoming 
waves,  if  the  Helmholtz  equation  had  been  obtained  by  making  the 
assumption 

?(Xf  /,  t )  -  1 

in  the  equation  for  acoustical  waves 

A» 

The  author  did  not  find  this  analogy  completely  convincing,  but 
we  shall  show  in  the  next  section  that  Eq.  (10)  can,  indeed, 
be  interpreted  as  a  superposition  of  outgoing  acoustical  waves. 

For  someone  willing  to  accept  this  interpretation,  there  is  no 
need  for  the  discussion  of  the  next  section. 

Next,  we  derive  solutions  corresponding  to  a  family  of 
particular  solutions  which  include  the  counterparts  of  Eqs.  (10) 
and  (11)  for  the  two  dimensional  case.  We  omit  the  z  dependence 
in  Eq.  (7)  and  introduce  polar  coordinates 
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0.)  One  then  obtains 


cos 

Jon.  f  m  6  ) 


This  is  Bessel's  equation.  One  then  obtains  the  solutions 

,  jj  .  f  COifnv6) 

l  6  / 


COS(nv  &  ) 
Jirof'**  &  ) 


The  expressions  H  behave  asympototically  as 


This  suggests  that  the  expressions  18  and  20  represent  outgoing 
and  the  expressions  19  and  21,  incoming  waves. 

Eq.  (15)  shows  that  g  is  oscillatory  for  r  >  m/y  and 
nonoscillatory  for  r  <  m/y.  (In  the  nonoscillatory  region,  g 
is  monotonic  except  for  possibly  one  minimum  of  the  absolute 


value  of  g.)  If  the  nonoscillatory  region  of  the  solution 
resembles  that  of  the  Laplace  equation,  the  region  increases  with 
m  and  decreases  with  increasing  y. 

We  derive  the  corresponding  particular  solutions  also  for 
the  three  dimensional  case.  Usually  one  needs  them  only  for  r'=0. 

We  introduce  a  system  of  spherical  coordinates  with  radius  r,  \p 
longitude  ,  and  0  latitude  and  set  £  =  sin  0 ,  where  the  z  axis  points 
in  the  north  direction;  that  is,  we  introduce 

r^=  x *■  ty*’  + 

u  /  (22) 
J*-  z/r 

arctf(y/x) 

then  one  obtains 

/xx  V**  T-<tr  +  £  tfY  (23) 


The  hypothesis 


fai.r)  -jC 


gives 


Next,  an  expression 

.  *<*  .  M>)  ,  .  /  - 

r  fr>f*.cV 

is  substituted  into  Eq.  (25)  .  Multiplying  the  result  by 
2 

r  /(f(r)g(c))  and  introducing  a  separation  constant  X,  one  then 


obtains 


jy  /"-SJ  ’ 


(21') 


Eq.  (27’)  is  the  equation  of  associated  Legendre  functions 
(spherical  harmonics)  usually  denoted  by  Pmn.  Notice  that 
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Eq.  (21')  does  not  contain  the  frequency  p.  The  separation 
constant  is  determined  from  Eq.  (21')  by  the  requirement  that 
the  function  finally  obtained  be  regular  at  the  poles  £  =  +1. 
One  obtains 

X  =  n(n  +  1) 

(This  explains  the  introduction  of  n  anticipated  in  Eq.  (26).) 
Moreover,  one  has 


m  <  n 


(see  for  instance  Reference  4,  Chapter  11) 
Equation  27  gives 

u  .fit.) 


cr.  —  , 

c/r*-  r  <fr 


n i 


Substitution  of  X  into 


(28) 


This  differential  equation  has  a  regular  singular  point  with 

exponential  n  and  -(n+1)  at  r  =  0  and  an  essential  singularity  at 

2 

infinity.  The  series  development  proceeds  in  powers  of  r  .  The 

/  f.  \ 

function  f '  ’  can  be  expressed  in  terms  of  elementary  transcendental 
functions  .  One  sets  for  this  purpose 


-fi***) 

r 


(29) 


then  one  obtains 


0  -  P0  -  *&£  - 

As. 

A  power  series  hypothesis  for  f 

A 


?-Z°k 

leads  to  the  recurrence  relation 

A'fl,  ~  k  ) 


hence 


h,  tv!  / 

fa-k)!  k! 


(30) 


(31) 


(32) 
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The  coefficients  are  zero  for  k  <  0,  and  for  k  >  n  for  the  factorials 
of  negative  numbers  are  infinite.  In  other  words,  the  function 

A 

f  is  a  polynomial.  A  second  solution  would  start  with  a  2n  +  1. 

It  is  then  given  by  an  infinite  series.  But  a  second  solution 
can  be  obtained  in  a  different  way  (as  we  shall  see) . 

For  a  steady  flow  (y  =  0)  ,  fn(r)  is  given  by  r~n_1  and  rn. 

In  the  far  field  only  the  negative  powers  play  a  role.  In  a  steady 

flow  it  is,  therefore,  possible  to  formulate  far  field  conditions 

(where  r  is  large)  in  terms  of  a  very  limited  number  of  such 

-1  -2 

particular  solutions,  usually  only  r  and  r  (poles  and  dipoles) . 
The  contributions  of  single  poles  vanish  if  one  has  closed  bodies. 

In  oscillatory  flows,  all  particular  solutions  decrease  for  large 
r  and  r  it  may  therefore  be  necessary  to  take  a  large  number  of 
these  particular  solutions  into  account. 

Eq.  (28)  shows  that  the  solutions  fn  are  nonoscillatory 
2 

for  (yr)  <  n (n+1) . 

The  differential  equation  for  fn  (Eq.  (28))  has  real 
coefficients.  It  is  therefore  satisfied  separately  by  the  real 
and  the  imaginary  parts  of  the  expression  Eq.  (29) . 

If  a^  =  1  is  real,  it  then  follows  from  Eq.  (30)  that 
is  imaginary.  This  means  that  the  development  of  the  imaginary 
part  of  the  expression  (29)  has  as  lowest  power  r  n.  However,  the 
exponents  of  the  regular  singular  point  r  =  0  are  -n-1  and  +n. 

It  follows  that  the  lowest  power  in  the  development  of  the  imaginary 
part  is  actually  r+n.  The  expression  (29)  contains,  accordingly, 
the  solution  pertaining  to  the  exponent  +n  as  well  as  that  to  the 
exponent  -n-1. 

For  n  =  0  one  obtains  the  fundamental  solution,  Eq.  (10) . 
Because  of  the  presence  of  the  exponential  function  in  Eq.  (29) , 
one  will  surmise  that  also  these  expressions  are  related  to  outgoing 
waves.  (Of  course  the  real  or  imaginary  parts  by  themselves  are  a 
mixture  of  incoming  and  outgoing  waves.) 
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Incoming  waves  are  obtained  by  taking  the  conjugate  complex 
of  the  expression  29.  Let  us  introduce  the  notation 

?  (33) 

4W>  „£(t)  *V**^/V'  > 

(34) 

'  r 

where  f  is  given  by  Eqs.  (29),  (30),  and  (31),  and  *  denotes  the 
conjugate  complex.  Solutions  with  last  subscript  2  represents 
incoming  waves.  One  notices  that  4>  (m,n,i,l)  £s  £jje  conjugate 
complex  of  n'  ^  ,  There  exist  orthogonality  relations 

-  / 

There  exist  also  formulae  for  the  normalization  constant,  but  they 
are  not  needed  here.  Besides,  one  has  orthogonality  relations 
between  the  trigonometric  functions  sin  m\|>  or  cos  m\J>.  Then,  one 
has  for  constant  r 

where  do  is  the  surface  element  of  the  sphere.  Notice  that  one 
may  have  k1  =  k2»  The  derivative  of  the  function  <j)(m»n»i»k)  w-th 
respect  to  r  is  the  derivative  with  respect  to  the  normal  of  this 
surface,  it  will  be  denoted  by  d^/dn.  Now  consider 

f  %  %  k)*_  ”'J  *'*  **  ^  ^  ^ 

This  expression  obviously  vanishes  for  (m^,  n^,  i^)  ^  (m2,  n2,  i2) 
If  these  subscripts  agree,  and  k^  and  k2  are  different,  k^  =  1 
and  k2  =2,  say,  then  one  obtains 

where  3x  is  the  spatial  angle  of  the  surface  element.  If  k^  =  k2 
1,  then  one  obtains  .. 

J I**'  K 
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The  expression  in  the  bracket  is  the  Wronskian  of  two  linearly 

dependent  solutions  of  Eq.  (28) .  Because  of  the  factor  2/r  of 

the  term  df  ^j'dr  in  this  equation,  the  Wronskian  is  given  by 
_2 

const  r  .  This  shows  that  the  integral  is  independent  of  r. 
Thus ,  one  has 

=  0J  c/-  «>,  t>,)  K  ) 

If  all  superscripts  are  the  same,  then  one  obtains  a  constant 
independent  of  r.  So  far  this  formula  holds  only  for  surfaces 
r  =  const.  In  Section  IV  it  will  be  extended  to  surfaces  of 
arbitrary  shape.  The  formula  can  be  used  to  decompose  a  given 
function  4>  into  incoming  and  outgoing  waves. 
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SECTION  III 

INTERPRETATION  OF  THE  SOLUTIONS  DERIVED  IN  SECTION  II 

We  start  with  the  equation  for  the  small  acoustical 
perturbations 

4>  +  <j>  +  4>  —  a  ^<j>  .  =0  ^ 

Yxx  Tyy  zz  rtt 

Let  r  and  r'  be  vectors  with  components  x,  y,  z  and  x' ,  y' ,  z', 
respectively,  and  let 

r  =  |r  -  r' |  (37) 

Then,  one  has  the  well  known  particular  solutions 

4>  (r ,  r',  t)  =  r  1(f(r-at))  (38) 

where  f  is  an  arbitrary  function.  Taking  for  f  a  function  which 
has  finite  support  (that  is  a  function  which  differs  from  zero  only 
in  a  finite  interval)  one  clearly  obtains  the  picture  of  a  wave 
package  which  travels  with  the  speed  a  in  the  direction  of 
increasing  values  of  r.  (The  shape  of  this  package  changes  because 
of  factor  r  ^.)  The  expression  (38)  can  be  interpreted  as  the  field 
caused  by  a  source  at  the  point  r  =  r',  its  strength  is  given  by 
the  momentary  mass  flux  through  a  small  sphere  around  this  point. 
Assuming  the  density  to  be  1,  one  obtains  the  source  strength  at 
time  t 

/(-at) 

The  total  mass  flow  during  a  time  interval  from  t^  to  t2  is  given 
by 

_*jT  J 

Choosing  for  f  a  delta  function,  one  obtains  the  effect  of  a 
momentary  source  at  (at  time  t') 

or 

fM  * 
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The  total  mass  flow  of  this  source  is  given  by 

O' ft t£'  /  SC*)  Sr  -  - 

tst~t  + 

The  potential  for  a  wave  caused  by  a  momentary  source  of  strength 
1  at  time  t  =  t'  and  at  a  point  r  =  r 1 ,  is  therefore  given  by 


-fify  Osir' £(r- a,ft-t')) 


(39) 


A  flow  field  which  can  be  represented  as  a  superposition  of 

such  waves  (with  different  times  t'  and  different  centers  r')  and 

where  r'  is  confined  to  a  finite  region,  satisfies  the  condition 

that  at  a  sufficient  distance  it  consists  only  of  outgoing  waves . 

This  is  our  basic  definition.  Now  consider  a  superposition  of  such 

waves  whose  strength  is  given  by  exp(ivt').  Then,  with  Eq.  (39), 

t 

<j>  (r,  *,+)  -  -fit)' 4-  r'f  M/a  (cut' )  «,ft-t‘))St' 

-  Co 

Let 

ym  r  -  Q-  ft  -  t' )  .t/v  *  ect/t 

The  integrand  is  zero  except  for  v  =  0,  because  of  the  presence  of 
the  delta  function,  that  is  for 


One  thus  obtains 

fifr]  t )  -  -frty'r'  tx/iftvft  (40)- 

This  result  can  be  obtained  directly  from  Eq.  (38). 

The  equation  for  an  unsteady  flow  linearized  for  the  vicinity 
of  a  parallel  flow  in  the  x  direction  with  velocity  U  arises  from 
Eq-  (36) ,  if  one  views  the  flow  field  from  a  coordinate 
system  traveling  with  velocity  U  in  the  negative  x  direction. 
Setting 

x  -x  +  Ut)  £  -  i-Zj  -t.t  (4i) 

y,  *  ”  f  (*>  £  J 

one  arrives  at  Eq. (1)  . 
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We  now  derive  particular  solutions  of  Eq.  (1)  by  a 
superposition  of  momentary  sources  of  strength  exp(ivt)  moving 
in  the  negative  x  direction  with  velocity  U  along  a  line  y  =  y1  = 
const,  z  =  z'  =  const  so  that  it  arrives  at  the  station  x  at  time 
t'  =0.  The  flow  field  so  arising  is  viewed  from  a  system  of 
coordinates  which  moves  with  the  same  speed  in  the  negative  x 
direction  so  that  the  origin  of  the  moving  system  of  coordinates 
arrives  at  the  origin  of  the  system  of  coordinates  fixed  in  the 
air  at  rest  t  time  0.  In  the  moving  system  the  source  lies  at 
x  =  x'  at  all  times.  Accordingly,  we  chose 


t'.O <'  -  Vi,  y,  *]* 


Let 


(43) 


Then  one  has 


'A 


.  /if.  r'/  .  [ff  7- *')  -  l/fi-t'j/  + 

Using  Eq.  (39),  one  then  obtains  for  the  potential  in  a 
moving  system  of  coordinates 

r’t )  m  f  ')  ' 


(44) 


(45) 


with  _  , 

^  -  f*j  &  *) 

-  >  /  r  /j  * 

r  '«  (x,  Vj*  J 

We  set 

rm  r-  4,(t- 1' )  -  ffrx-x')- ItCi-t'j)*’-*  £  / 


(46) 


c/v  «  fr  Uft-t  )JV+  cvje/t 

Hence  t 

Aff,  r;t).-a  (**/' J Kvi')[Vh-*') -  V'ft-ff* “ rj  ?(»</» 


(47) 


The  integrand  is  zero  except  for  v  =  0.  Hence,  from  Eq.  (46) 

t-f- 


(48) 


Only  the  positive  root  appears  because  of  the  limits  of  the  integral, 
Eq.  (45) .  For  v  =  0  one  obtains  from  Eq.  (46) 

r  =  a  (t  -  f) 


One  obtains  for  the  term  in  the  bracket  of  Eq.  (47) 

JJ(x-k)  -  v*(t-t')+  o_r  ~  [(, x-k')x+  //-  1L)$ V  /x 

Using  Eq.  (48)  ,  we  finally  obtain 

x  [(t-X^f 

This  expression  is  singular  at  the  point  r  =  r ' .  It  leads  directly 
to  the  particular  solution,  Eq.  (10).  This  shows  that  this 
particular  solution  represents  the  potential  of  an  oscillating 
source  moving  in  air  at  rest.  Thus,  the  fundamental  solution, 

Eq.  (10) ,  represents  a  superposition  of  outgoing  acoustical 
waves . 


The  corresponding  solutions,  Eq.  (18)  ,  can  be  obtained  by 
the  so  called  method  of  descent;  tnat  is,  one  considers  a  distribu¬ 
tion  of  sources  of  the  form,  Eq.  (10)  ,  which  is  constant  along  a 
line  x  =  x',  y  =  y'  of  the  three  dimensional  space.  The  potential 
is  then  obtained  by  an  integration  with  respect  to  z"' ,  the  result 
is  independent  of  z.  One  obtains  in  this  manner  one  of  the  integral 
representations  of  the  Hankel  function. 

Differentiating  the  three  dimensional  fundamental  solution 
Eq.  (10),  with  respect  to  either  x',  y'  or  z'  or  equivalently 
with  respect  to  x,  y,  or  z,  n  times  one  obtains  expressions  of  the 

A 

form,  Eq.  (29)  ,  in  which  the  functions  f  are  polynomials  of  the 
degree  n.  Such  expressions  satisfy  the  original  Helmholtz 
Eq.  (7)  . 
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Not  all  of  these  derivatives  are  independent  because  the 
function  to  be  differentiated  satisfies  the  differential  Ea.  (7). 

The  process  of  differentiation  is  the  limit  rf  a  superposition 
of  expressions  (10)  .  Accordingly,  these  derivatives  amount  to  a 
superposition  of  outgoing  acoustical  waves. 

If  one  carries  out  a  fixed  number  of  n^  of  such  differentiations, 
then  one  obtains  a  finite  number  of  linearly  independent  expressions 
of  the  form 

j>  *  f(r> 

where  the  function  f(r)  has  the  form 


ffy)  -=  -I*fs(-Lpr)  j>(Y) 


-  *%,.  -/ 


in  which  the  function  P(r)  is  a  polynomial  of  degree  n^. 

These  particular  solutions  can  be  represented  by  the  4>m,n'1,k 
defined  in  Eq.  (33).  At  infinity  the  solutions  with  exp(-iyr) 
and  exp(+iyr)  (see  the  definition  of  f  in  Eq.(39)  )  behave  as 

r  1.  Therefore,  one  must  admit  functions  <j)m'n,i'k  with  superscript 
k  =  2  as  well  as  k  =  1.  Since  at  the  origin  one  of  these 
derivatives  has  as  lowest  powers  r-nl-1f  one  needs  to  admit 


particular  solutions  «j>  with  m  <  n  <  n^. 


There  are  only  a  finite 
number  of  such  solutions.  Then  one  has  the  following  situation. 
Consider  ip  and  0  as  fixed.  Along  such  a  ray  through  the  origin 
the  function  <f>  obtained  by  forming  the  mixed  derivative 

has  the  form 

....  /  ffy?) 


r  ' '  £ 

The  particular  solutions  (Eq  .  (33))  used  for  the  representation 
of  this  function  combine  into  an  expression  of  the  form 

r’*~> +  rK'  ' 

Equating  these  two  expressions,  one  obtains 


r*  * 7  /  r ;  -  r’*  fry)  +  r  v  ^ 


and  hence 


-Hf-/  ~ft) 


Ve/ur  ) 


r"  -  2 


2*}  ) 
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Now  one  has  a  polynomial  in  r  on  the  left  and  a  transcendental 
function  of  r  on  the  right.  This  equation  can  be  satisfied  only 

/  3  \ 

if  P  ni  =0*  In  other  words,  the  representation  of  the  derivatives 
contains  only  functions  4>  ^m,n'  ^  in  which  the  last  subscript  is 

k  =  1.  Now  it  is  a  matter  of  counting  the  number  of  linearly 
independent  derivatives  and  the  number  of  functions  4> (m»n» i • I) 
order  to  demonstrate  that  the  relations  so  constructed  can  be 
inversed;  in  other  words,  that  the  functions  ^ (m » n # i # 1 )  can  be 
represented  as  a  linear  combination  of  derivatives  of  the 
fundamental  solution  (Eq.  (10)).  Then  it  follows  that  the 
expressions  (Eq.  (33))  with  last  subscript  1  can  be  interpreted 
as  a  superposition  of  outgoing  acoustic  waves. 

Specific  formula  are  not  needed  in  the  present  context. 

They  are  derived  in  the  appendix. 


Now  one  has  a  polynomial  in  r  on  the  left  and  a  transcendental 
function  of  r  on  the  right.  This  equation  can  be  satisfied  only 
if  P^3ni  =0*  In  other  words,  the  representation  of  the  derivatives 
contains  only  functions  <p  ^m'n'  in  which  the  last  subscript  is 

k  =  1.  Now  it  is  a  matter  of  counting  the  number  of  linearly 
independent  derivatives  and  the  number  of  functions  4>  ^ra,n'  1 '  in 
order  to  demonstrate  that  the  relations  so  constructed  can  be 
inversed;  in  other  words,  that  the  functions  <J>  ^m,n' can  be 
represented  as  a  linear  combination  of  derivatives  of  the 
fundamental  solution  (Eq.  (10) ) .  Then  it  follows  that  the 
expressions  (Eq.  (33))  with  last  subscript  1  can  be  interpreted 
as  a  superposition  of  outgoing  acoustic  waves. 

Specific  formula  are  not  needed  in  the  present  context. 

They  are  derived  in  the  appendix. 
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SECTION  IV 
GREEN'S  THEOREM 


So  far  we  have  derived  a  system  of  particular  solutions  of 
the  Helmholtz  equation  including  a  fundamental  solution  and  we  have 
shown  that  they  can  be  interpreted  as  a  superposition  of  outgoing 
acoustical  waves.  These  particular  solutions  satisfy  the  far  field 
conditions.  Flow  fields  which  can  be  represented  as  a  superposition 
of  such  expressions  will  also  satisfy  the  far  field  conditions. 

This  statement  can  be  regarded  as  a  characterization  of  such  flow 
fields.  In  a  further  development,  we  bring  this  criterion  into 
a  form  which  is  more  practical.  This  is  done  by  means  of  Green's 
formula . 

Let 

L(i)  +<}yy  (49) 

and  let  tp  and  oj  satisfy 


/  (j>)  -  hfr  )  -  0 

i(u>)  =  0 


(50) 

(51) 


where  h(r)  is  considered  as  known.  Let  R  be  some  region  in  the 
x,y,z  space  and  &R  its  boundary.  Within  this  region  Eqs.  (50) 
and  (51)  are  satisfied  everywhere.  Then  one  has,  because  of 
Eq.  (51) 

///  [  L  (<f>(r‘J)  -  h(r')J  tof  r'jjtf'  =  o 
ft 

where  dv'  is  the  volume  element  in  the  three  dimensional  space 
of  the  variable  x',y',z'.  Denote  by  ipn,  and  the  derivative 
in  the  direction  of  the  outer  normal  in  the  x',  y',  z'  space. 
Then  one  has  the  familiar  formula 

-  jf// ft r?)  *  * 

Taking  for  w(r')  the  expression  o)(f,?')  given  in  Eq.  (10) 
surrounding  the  singular  point  r  =  r'  by  a  small  sphere  in  the 
r'  space,  one  obtains  in  a  familiar  manner 
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PT 


f(r)  -  III-  +  UJ?  r‘)  fifr')]  </o'  +f/fh(r')U>{r  9')Jy '  (52) 

dR  * 

For  our  application,  one  considers  a  region  bounded  by  the  two 
surfaces  91^  and  9r2 .  The  normal  derivatives  always  refers  to  the 
outer  normal  of  the  region  R.  One  is  not  justified  to  disregard 
the  contribution  of  the  surface  91*2  even  if  it  should  move  to 
infinity,  for  the  expressions  ip  and  u>  do  not  die  out  sufficiently 
fast.  The  limiting  process  is  justified  for  the  steady  case 
(u  =  0) .  Eq.  (52)  is  an  identity.  It  holds  only  if  the 
values  of  <t>n»  (r')  and  <p(r')  are  taken  from  a  function  <J>  which 
satisfies  Eq. 
expression 


(50) .  For  a  further  discussion,  consider  an 


(53) 


The  function  u)(r,r')  satisfies 


L  r  )  =  V 

where  the  subscript  r  under  L  indicates  that  the  operation  L  is 
to  be  carried  out  with  respect  to  the  variable  r,  while  r'  is  a 
fixed  parameter.  Because  of  the  normalization  of  <i)(r,r'),  one 
then  finds  from  the  last  formula  that 


lljft!)  *  hfr> 
everywhere  inside  and  R  and 
L(<f>Cr))  =  O 


(54) 


(55) 


outside  the  region.  For  small  r  =  |r  -  r'|  the  function  u(r,r') 
behaves  as  -4Tr/r;  that  is,  exactly  as  the  fundamental  solution 
of  the  Laplace  equation.  The  surface  integrals,  therefore 
represent  a  layer  of  poles  and  dipoles  with  densities  given 
respectively  by  f^fr')  and  f2(r').  A  layer  of  poles  generates  a 
function  <J>  which  is  continuous  everywhere  and  for  which  the  normal 
derivative  has  a  jump  given  by  the  density  of  the  poles,  here  by  f^. 
A  layer  of  dipoles  oriented  in  the  direction  of  the  normal  gives 
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a  function  tf>  for  which  the  normal  derivative  is  continuous  everywhere 
and  for  which  the  jump  of  the  potential  is  the  negative  of  the 
dipole  density  (-f2)  •  This  is  an  interpretation  of  Eq.  (53) . 

If  f^  happens  to  be  equal  to  the  values  of  <J>n  pertaining  to  an 
actual  solution  and  -f2  happens  to  be  the  value  of  4> ,  then  the 
potential  outside  of  R  is  identically  equal  to  zero  and  one  recovers 
inside  the  function  <Mr)  . 

Let  9R  be  the  contour  of  a  body.  An  expression 

f(?)  .  // 1-*?,?), f,(ri  -  o'  ,  ff/kt?') 

it  f 

satisfies 

Lfy)*  hfr) 

and  the  far  field  condition  for  <J>,  for  it  is  a  linear  combination 
of  functions  which  satisfy  this  condition.  If  f2(r')  and  f^(r') 
happen  to  be  identical  with  <j>  and  its  normal  derivative  at  9R, 
then  the  above  formula  gives  a  representation  for  4>  in  the  field. 
(Usually  <J>  is  prescribed  and  4>  is  unknown.)  Notice  that  in  this 
derivation,  we  have  used  results  of  potential  theory  rather  than 
Green's  formula.  A  comparison  with  Eq.  (52)  shows  that  the 
integral  over  the  outer  boundary  does  not  appear.  How  this  comes 
about  is  seen  by  the  following  argument. 

Let  h(r')  =  0.  Then  one  has 

JJ  /-  / ,  fr)  *>>(?')  +</>  fr'j  tvjfrjje/o  '  +  //[  fa f*J  u/?')Jdo ' 

This  shows  that  the  integral 

// [-fa  y  f  (*’)%'»  fr'tj Jo' 

i/t 

remains  unchanged  if  one  carries  out  a  continuous  deformation  of 
the  surface  9R  through  a  region  in  which 

L  ( (p  )  =  o  tnd  o 

Notice  that  (f>  or  w  need  not  satisfy  special  far  field  conditions. 

Such  a  deformation  can,  for  instance,  be  carried  out  for  the 
spherical  surface  over  which  the  integration  in  Eq.  (34)  is  carried 
out.  We  apply  this  result  to  two  functions  co(r,r')  and  w(r",r') 
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lie  inside  of  the 


Figure  2.  Surface  9R  in  the  r'  -  space. 

The  surface  9R  can  be  contracted  into  two  small  spheres  around  the 
singular  points  r'  =  r  and  r'  =  r",  and  then  the  integrals  can  be 
evaluated.  One  finds  that  the  contributions  of  the  two  spheres 
cancel  each  other.  Thus, 

JJ  [  (otf  rj  rj  -  v,  fr,  r'j  coft*  r'j]  </o',  e  (56) 

A  function  <J>(r)  that  satisfies  the  far  field  condition  can  be 
represented  as 

<f>(r)  ~  %  fig,  %*) 

This  function  evaluated  at  r  =  r'  gives 
tfr')*  £  fa  vfr;?/) 

Thus,  one  finds  that 

f/  [ =o  (57) 

This  explains  why  the  integral  over  the  outer  contour  is  not  present 
in  Eq.  (52),  if  the  far  field  conditions  are  satisfied.  It  is, 
of  course,  assumed  that  h(r')  =  0  outside  of  the  outer  contour, 
for  Eq.  (57)  does  not  hold  if  one  of  the  singular  points 
r'  =  r  or  r'  =  r"  lies  outside  of  the  surface  of  integration. 

The  relation  (Eq.  (57))  constitutes  a  far  field  criterion. 

It  must  be  satisfied  for  every  choice  of  r  (although  one  will  then 
obtain  redundant  conditions) .  It  therefore  holds  also  for  derivatives 
of  this  equation  with  respect  to  r,  that  is,  if  the  function  u)(r,r') 
is  replaced  by  its  derivatives  with  respect  to  x,  y,  or  z, 
usually  at  the  origin.  But  these  derivatives  can  be  expressed 
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by  the  functions  <}> 
a  further  relation 

<("»,  -nj  ijt)M 


(n,  m,  i,  1) 


(Eq.  (33)).  Thus,  one  obtains 
*•>  fJM 


for  all  functions  <f> '  ’  '  '  ;  and  an  arbitrary  outer  boundary 
surface  of  the  region  under  consideration.  One  can  gain  from 
Eq .  (57)  another  form  of  the  far  field  criterion.  Let  S  be  a 
surface  in  the  flow  field  which  includes  all  points  at  which 
h(r')  f  0  but  lies  inside  of  the  boundary  3R  of  the  computed 
field.  Let  <j>  be  the  solution  of  a  problem 


(58) 


L(<f>  (r))  -  hfr)  -  c 

which  satisfies  the  far  field  conditions  and  let  <j>  be  the  values 

s 

of  <f>  at  the  surface  S.  Then  one  can  solve  the  Dirichlet  problem, 

which  has  L<f>  =  0  inside  of  S  (including  the  region  within  an 

inner  bounding  body  surface) .  This  can  always  be  done  provided 
2 

that  y  is  not  an  eigenvalue.  This  means  that  L4>  =  0  with  the 
boundary  condition  <f>  =  0  at  S  does  not  have  a  nontrivial  solution. 
Excluding  this  possibility,  one  has  now  a  function  <(>  which  is 
represented  outside  of  S  by  the  solution  <{>  of  the  original  problem 
satisfying  the  far  field  condition  and  inside  of  S  by  the  solution 
of  the  Dirichlet  problem  just  described.  At  S,  this  function  <p  is 
continuous,  but  it  has  a  jump  of  the  normal  derivative.  This  flow 
field  can  then  be  represented  by  a  layer  of  poles  distributed  over  S 
where  the  pole  density  is  given  by  the  jump  of  the  normal  derivative. 

Thus,  we  have,  as  the  representation  for  the  original  <j>  outside  of  S 


<f>fr)  -  Jf  f-fr')  u> (■?,?')  e/ o'  ,  r  c  S  (59) 

S 


This  representation  can  be  applied  in  particular  to  the  outer  surface 
3R  where  the  farfield  conditions  are  to  be  formulated.  Here  the 
function  4>  is  represented  by  a  superposition  of  particular 
solutions  cj(r,r')  which  satisfy  the  far  field  conditions  and  where 
the  singular  points  are  arranged  at  a  surface  S. 
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The  limiting  case  where  the  surface  S  approaches  the  outer 

boundary  of  the  computed  flow  field  is  admissible  (and  sometimes 

2 

practical) .  In  cases  where  y  is  an  eigenvalue  of  the  problem 
for  the  function  S,  or  even  close  to  an  eigenvalue,  one  will  find 
that  the  function  f  in  Eq.  (59)  becomes  very  large.  For  numerical 
reasons,  such  cases  must  be  excluded. 


SECTION  V 

FAR  FIELD  CRITERIA 

We  started  from  the  definition  that  any  superposition  of 
expressions  which  ultimately  can  be  considered  as  a  linear 
combination  of  expression  representing  outgoing  acoustical  waves 
satisfies  the  far  field  conditions.  From  this  idea  the  following 
criteria  have  been  derived. 

1.  The  solution  is  a  superposition  of  expressions  (33) 

2.  The  solution  satisfies  the  condition  Eq.  (58) 

3.  The  solutions  can  be  represented  in  the  form  of 
Eq .  (59) 

4.  The  Klunker -Traci  formulation 

Formulations  1  and  3  are  closely  related  to  each  other.  In  each 
case  the  solution  at  the  outer  boundary  of  the  computed  flow 
field,  (this  includes  as  well  as  <J>)  is  represented  by 
a  linear  combination  of  particular  solutions  which  satisfy  the 
far  field  conditions. 

In  the  formulation  1,  there  arises  the  question  of  convergence 
if  the  outer  surface  of  the  region  under  consideration  deviates 
strongly  from  a  circle. 

2 

In  the  formulation  3,  one  has  the  restriction  that  y  must 
not  be  an  eigenvalue  of  the  Dirichlet  problem  for  the  chosen 
surface  S.  We  shall  see  in  an  example  that  this  may  indeed  be 
an  impediment.  In  the  formulation  2,  the  convergence  problem 
appears  in  a  different  form.  The  functions  <d  play  the  role  of 
a  test  function.  If  the  value  of  r  varies  rather  strongly  along 
the  outer  surface,  and  if  n  is  large,  then  because  of  the  character 
of  the  function  f^n^(r),  the  functions  are  large  for  small 
values  of  r  and,  comparatively  small  for  large  values  of  r.  For 
large  values  of  n,  the  particular  solutions  behave  roughly  as 
r~n.  This  means  that  the  particular  solutions  <f>  become  nearly 
linearly  dependent  for  large  values  of  n  ,  if  r  varies  strongly 
along  the  outer  surface. 
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We  snail  show  in  an  example  that  the  exclusion  of  the 
vicinity  of  eigenvalues  in  the  application  of  the  criterion  3 
may  be  necessary.  If  y  is  an  eigenvalue,  then  there  exists  one 
or  more  functions  <p  for  which  L$  =  0  in  the  interior  of  S  and  for 
which  <p  =  0  at  the  surface  S.  The  normal  derivative  of  this 
function  in  the  interior  is,  of  course,  not  zero.  This  function 
can  be  continued  outside  of  S  by  4>  =  0.  A  layer  of  single  poles 
with  density  given  by  the  values  of  $n  pertaining  to  the  eigen¬ 
function  will  represent  this  function  4>  and  therefore  give  $  =  0 
outside.  In  other  words,  one  has  one  or  more  density  distributions 

which  will  not  contribute  to  the  functions  d>  and  $  at  the  outer 

n 

edge  of  the  computed  flow  field.  We  shall  see  in  an  example  that 
the  particular  solutions  for  4>  which  are  lost  because  of  this 
phenomenon  are  needed. 

The  formulation  2  of  the  far  field  conditions  is  particularly 
useful  if  one  has  to  make  a  refinement  to  an  existing  approximate 
solution  and  one  knows  that  the  contributions  of  long  waves  are 
small.  (These  are  the  solutions  for  which  one  may  find  oneself  in 
the  vicinity  of  an  eigenvalue.)  The  criterion  2  works  rather  well 
for  long  waves  (low  values  of  n)  and  the  criterion  3  for  short  waves 
(large  values  of  n) .  In  this  sense  the  criteria  2  and  3  are 
complementary.  Further  remarks  about  the  application  of  this 
criterion  will  be  made  later. 

The  formulation  of  Klunker  and  Traci  is  based  on  Eq.  (53) . 

Here  f^  is  the  given  normal  component  of  g>  and  f2  the  unknown 
potential.  In  essence,  this  criterion  is  rather  similar  to 
criterion  3  shown  above  except  that  the  surface  S  is  contracted 
to  the  body  surface  and  that  as  a  consequence,  one  takes  the 
presence  of  sources  h(r')  within  the  field  into  account. 

There  is,  however,  a  difference.  In  the  formulation  3,  the 
function  f  is  considered  as  unknown  and  one  expresses  4>  and  <j>n 
at  the  outer  edge  of  the  flow  field  in  terms  of  this  function. 

In  the  final  computed  flow  field  it  must  be  possible  to  represent 
<J>  as  well  as  <J>  by  means  of  a  suitably  chosen  function  f.  The 
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same  can,  of  course,  be  done  with  the  use  of  Klunker's  formulation. 
In  this  case  the  potential  at  the  surface  plays  the  role  of  the 
function  f,  but  one  does  not  make  direct  use  of  the  fact  that  this 
function  has  this  specific  interpretation. 

In  the  Klunker JTraci  formulation  one  proceeds,  howeve^ 
in  a  different  manner.  One  uses  the  expression  (Eq.  (52))  to 
compute  only  $  at  the  outer  boundary.  Next,  one  uses  these  values 
of  4>  as  the  far  field  condition  and  computes  the  flow  field.  In 
this  computation,  one  finds  the  surface  potential  <j> .  One  has  the 
desired  solution  if  the  surface  potential  so  found  agrees  with 
the  potential  used  to  compute  4>  at  the  far  field  boundary. 

The  difference  between  the  formulation  3  and  4  lies  in  the 
fact  that  in  formulation  4  one  disregards  the  condition  for 
at  the  far  boundary  which  one  could  obtain  from  Eq.  (52) ,  and 
replaces  it  by  the  requirement  that  the  function  f 2  which  appears 
in  Eq.  (53)  is  identical  with  the  surface  potential  obtained 
from  the  flow  field  computation.  The  Klunker  Traci  formulation 
is  well  suited  to  an  iterative  approach  in  which  one  updates 
the  surface  potential  in  each  iteration  step.  It  becomes 
inconvenient  if  iterations  should  fail. 
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SECTION  VI 
AN  EXAMPLE 


The  following  example  shows  how  these  conditions  can  be 
applied.  It  is  extremely  simple;  in  fact,  all  steps  which 
ordinarily  require  numerical  techniques  are  carried  out  analytically. 
The  sole  purpose  is  to  show  which  steps  must  be  taken  if  one 
applies  different  criteria.  We  shall  consider  solutions  of  the 
Helmholtz  equation,  interpreted  as  a  description  of  the  propagation 
of  small  perturbations  in  air  at  rest.  The  boundary  of  the  body  is 
given  by  a  sphere  which  pulsates  harmonically  in  time.  The  three 
dimensional  case  has  been  chosen  because  of  the  simplicity  of  the 
fundamental  solution.  The  field  has  spherical  symmetry  because  the 
boundary  conditions  are  assumed  to  have  spherical  symmetry. 
Accordingly,  we  deal  with  a  one  dimensional  problem.  It  is,  of 
course,  important  that  the  steps  for  which  we  use  analytical 
formula  must^in  a  realistic  case, be  carried  out  by  numerical 
techniques  in  one  or  three  dimensions.  While  in  the  present  case, 
the  far  field  conditions  depend  upon  only  one  parameter,  one  will 
have  in  reality  a  great  number  of  such  parameters. 

The  field  for  which  the  computations  are  carried  out  lies 
between  the  radius  Tq  (radius  of  the  pulsating  sphere)  and  a  radius 
r^,  at  which  the  far  field  conditions  are  to  be  applied.  At 
r  =  rQ  the  value  of  $  =  <j>rQ  =  const  is  prescribed.  The  differential 

equation  to  be  satisfied  reduces  to 

+  *  rfr  *  °  (60) 

The  fundamental  solution  which  satisfies  the  radiation  condition 
is  given  by 

p-  -  ft*)'*'  (61) 


The  problem  has  the  following  solution 

-5  K 


/  ■+  ro 


y.  [-(jf'fr-Tc)] 


(62) 
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It  is  readily  verified  that  it  satisfies  the  boundary  condition 


Because  of  the  spherical  symmetry  of  the  problem,  the  density 
functions  for  the  distribution  of  poles  or  dipoles  occurring  in 
various  formulae  will  also  be  spherically  symmetric.  At  various 
occasions,  for  instance  in  the  expression  (50)  we  shall  encounter 
expressions  4>(r)  =  // F(r')  (r-r ' ) do'  ,  where  the  integration  is  to 

be  extended  over  a  sphere  r'  =  const  and  where  F(r')  =  const.  Thus, 
one  has  to  evaluate  integrals 


<J> '  (r)  =  //  oj (r-r  '  )do  ' 


The  numerical  evaluation  which  is  usually  needed  is  replaced  here 

by  an  analytical  approach  in  which  we  construct  functions  which 

have  the  same  properties  as  these  integrals.  To  be  specific,  the 

functions  which  arise  by  these  integrals  satisfy  the  Helmholtz 

equation  inside  and  outside  of  the  spherical  layer.  At  infinity 

they  satisfy  the  radiation  condition  and  in  case  of  a  layer  of 

single  poles  with  density  1,  one  has  continuity  of  the  potential 

and  a  jump  of  the  derivative  normal  to  the  layer.  In  the  case  of 

a  dipole  distribution,  one  has  continuity  of  the  normal  derivative 

and  a  jump  of  -1  of  the  potential.  Let  the  potential  of  a  single 

layer  at  radius  r  be  given  by  <t>^,  and  the  potential  of  a  layer 

(2? 

of  dipoles  by  <J>  .  The  direction  of  the  outer  normal  which  must 

be  taken  at  r  =  rQ  is  the  negative  r  direction.  Then  one  has 


j>"‘~  -Jr  2 

tfr )*>«,£. r. 

One  verifies  that 


+  o)  -  -  o  )  -  O 


r  <r0 

r  > 


(63) 
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Moreover, 


fi*>u  ~  -^r  y~  j 

4>(*W )  =  f  ro)  -  lr  **/>{- j 

One  verifies  that 


y  <r0 

r>rc 


(64) 


/fy+o)  -6)  •  ' 

<j>“(r0+t)  -tZfa-o)  -  /“fr-*)  ~  /rfo+oJ  -  0 


Also  needed  in  the  analysis  is  a  general  solution  which  satisfies 
the  boundary  conditions  at  the  surface  of  the  body.  In  a  numerical 
approach,  one  usually  does  not  encounter  such  an  expression  directly 
The  analytical  formulation  is  used  to  replace  information 
usually  obtained  by  integrating  the  flow  equations.  One  has 


ro  (fro 


/-tl 


ir° 


-JUty>(-£ffr-rt)  +  c  y  jA.ro  (65) 


The  first  part  is  the  analytic  solution  of  the  problem  (which 
satisfies  the  inhomogeneous  boundary  conditions  at  r  =  r^,  the 
second  part  has  an  r  derivative  0  at  r  =  r^. 


After  these  preparations,  we  can  apply  different  far  field 
conditions.  In  the  formulation  one  expresses  the  far  field 
(that  is  the  field  at  r  =  r^)  by  a  superposition  of  solutions  (33) . 
Because  of  the  spherical  symmetry  of  the  present  problem,  only 
one  of  the  solutions,  Eq.  (33)  namely, 

')  m  -  {bxf' r' txpf-tfr) 


is  encountered,  with  a  coefficient  c^,  say. 
Then  one  has  at  r  =  r.^ 


<j>  (rf  )  *  -  Cf  f+tf'  Tr~'  if-  r,  ) 
4>Jy,)  -  CfWJ’Y  ^ 


(66  ) 


The  expression  Eq.  (65)  gives 

u  * 


3 


E 


*  C[~  -r,))  +j*r0  tos^{rrrv))]  +  ^^[cei^rrrc)) 

-f.r^fon,  -  rj)]  f 


Usually,  one  avoids  the  construction  of  a  general  solution 
by  numerical  solution  of  the  difference  equations.  (This  would 
require  far  too  many  free  parameters.)  In  the  present  problem, 
in  difference  form,  one  obtains  a  linear  system  for  the  values  of 
4>  at  the  points  of  a  grid.  If  one  eliminates  the  unknowns  from 
this  linear  system,  starting  at  the  profile  (here  r  =  rQ)  and 
proceeding  toward  the  outer  boundary  one  finally  obtains  relations 
between  the  values  of  tf  at  the  outer  boundary  and  at  the  row  of 
grid  points  next  to  it.  These  relations  are  equivalent  to  a  linear 
relation  between  <j>  and  (j>  at  the  outer  boundary.  In  the  present 
case,  the  corresponding  information  would  be  obtained  by  eliminating 
C  from  Eqs .  (67)  and  (68).  Into  this  equation,  one  then 

substitutes  the  expressions  Eq.  (66)  which  express  the  far  field 
conditions.  One  thus  obtains  an  equation  for  C-^.  The  solution 
can  obviously  be  found  by  inspection.  It  is  given  by  Eqs.  (67)  and 
(68)  with  C  =  0. 


In  the  second  formulation,  one  uses  Eq.  (58)  because  of  the 
spherical  symmetry  one  needs  to  consider  only  the  one  function 
<J>  (1 ' 1 ' 1 ' 1)  .  This  gives  the  relation 


-f-  tf**1/)  frr/  )  ~  0  (69) 


This  is  obviously  equivalent  with  Eqs.  (66)  .  This  equation  is 

now  combined  with  the  relation  between  d>  and  1  obtained  from  the 

n 

flow  field  (Eqs.  (67)  and  (68)  with  C  eliminated) .  The  problem 
is,  of  course,  very  similar  to  the  previous  case,  except  that 
the  constant  does  not  occur.  One  has  in  Eq.  (69)  a  direct 
relation  between  the  values  of  <J>  and  <f>n  at  the  outer  contour. 
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The  third  form  of  the  criterion  gives,  for  the  present 
example,  the  same  procedure  as  the  first  one  since  there  is  only 
one  function,  u)(r,r')  which  is  spherically  symmetric  and  satisfies 
the  far  field  conditions. 

We  describe  the  fourth  criterion  (Klunker-Traci)  in  a  manner 
which  is  suited  for  a  noniterative  approach.  First,  one  expresses, 
by  means  of  Green's  formula,  the  values  of  <|>  and  <j>r  at  r  =  r^  in 
terms  of  the  given  function  and  of  the  unknown  surface 


potential  <J>(rQ)  here  denoted  by  f^.  Making  a  token  evaluation  of 
Eq.  (52)  (actually  using  Eqs.  (63)  and  (64))  and  remembering 


that  at  r  =  r 


=  -<(> 


*0‘ 


one  obtains 


f>(r,)  *  v  +  f,[ ro )  -  (-  J 


^r0 


(70) 


frtr,  fW, 


1) 


Taking  the  original  form  of  the  Klunker  Traci  formulation, 
we  now  construct  a  solution  (in  practice  by  solving  the  difference 
equation  for  the  flow  field)  which  assumes  the  values  of  <J>(r-^)  at 
the  outer  edge  of  the  flow  field. 


still  contains  the  unknown  function  f^.) 


(Notice  that  this  formulation 
This  is  done  here  by 
choosing  the  constant  C  in  the  expression  (67)  .  It  is  useful  to 
rewrite  the  last  equation.  We  know  that  the  ultimate  solution  is 
obtained  for  C  =  0.  The  function  f^  pertaining  to  it  is  then 
given  by 


Now  we  set 


fi- 


'  fro 


t-U. 


pr0 


^ e>  (Pro 


(72) 


Then  one  obtains  from  the  boundary  condition  (70) 
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Hence, 


(  73) 


The  first  term  in  this  equation  agrees,  of  course,  with  Eq.  (62) 
and  the  first  term  in  Eq.  (65) .  To  satisfy  the  boundary 
conditions,  one  now  obtains  the  condition  (from  Eq.  (65)  and 
the  last  equation) . 


S'f,  [cot  jpn)-  ?  )  “  C  £  +fr.C*fc,(rrr.))] 

Hence, 

c  -  <r/  A  r° C0i^ r<>/>  ~  r*  ^ r>  > 

'  fr0[^^{rrr0))+frbC*f^rro))]  (  74  > 


Inserting  this  value  of  C  into  Eq.  (65),  one  obtains  an  expression 
for  <(>  which  satisfies  the  boundary  condition  at  r  =  r^  and  assumes 
the  values  of  <f>  at  r  =  r^  given  by  Eq.  (7  0)  .  (The  derivation 
shown  here  replaces  a  numerical  integration  for  the  flow  field.) 
This  solution  still  depends  upon  the  unknown  6f^.  From  this 
solution,  one  obtains  the  potential  at  the  surface  r  =  rQ. 

One  has  from  Eq.  (65) 

Afc)  -  -r°^o  ,  [>r  A  r<>  r* } ~  .  j  ^ 

/v^(70  <?-  ro  [toovfrl r,--r6)>  >1  cot{f.( rrr0)) 

Now  one  has  the  requirement  that  the  function  f^  originally 
assumed  agree  with  the  value  of  the  potential  at  r  =  r^  found 
from  integrating  the  flow  equation.  Hence, 

-  fr  jro  fr  -  r*4ro  y  ff  rc  )J,K/>  r,  ) 

'  JU,fc(rrrj)  +fr0c*/f.Cr,-rh)) 

Hence 

sir,.  i  .  0 

r>  -yi))  -r0)  ) 
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This  is  solved  by  6f^  =  0  (as  expected)  unless  the  factor  of  6f^ 
vanishes.  In  this  factor  some  simplifications  are  possible. 

One  ultimately  obtains 


<r/ 

.  ■  >  /  ..  /  / 


re  ) 


(75) 


+  fr0  Ctn^(r,-r0)) 

The  procedure  fails  if  sin  (yr^)  =  0.  Then  the  condition  is 
satisfied  for  any  value  of  6f^. 

In  an  iterative  procedure,  one  starts  with  an  assumed 
value  of  cj)  at  the  surface  r  =  r^,  that  is  with  some  choice  of 
6f^  and  then  computes  by  means  of  Green's  formula  the  potential 
at  r  =  r^.  According  to  Eq.  (73)  ,  it  deviates  crom  the 
unknown  potential  of  the  exact  solution  by 

With  this  boundary  condition,  one  integrates  the  flow  equation, 
and  arrives  at  a  surface  potential  whose  deviations  from  the 
exact  solution  is  given  by 


cai^>r0)  -  si+»fprJJ 
1  JvH/U*/r,  ~r-e))  -4  <?'rc  -*•)) 


-t  <y'rc  “H/f'fr, 

The  potential  so  obtained  is  then  taken  as  the  starting  point 
of  a  new  iteration  step.  The  iterations  converge  (in  this 
particular  case)  if  the  absolute  value  of  <5f^  decreases,  that  is, 
if 

I  £urc  Ct>o^MreJ  -  tv»(g*rc)]  j 

I  r„  -n))  I 


<  / 


In  general,  this  condition  is  not  satisfied.  It  is,  however, 

satisfied  for  y  sufficiently  small.  For  y  =  0  the  correct  value 

of  <p  is  obtained  in  one  step.  The  denominator  in  the  last 

2 

expression  vanishes  if  y  is  an  eigenvalue  of  the  problem 


fix  *  fa  * 
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with  boundary  conditions  <t>  =  0  at  the  surface  of  the  body  and 

<P  =  0  at  the  outer  contour.  If  one  is  in  the  vicinity  of  such  an 
eigenvalue,  and  one  tries  to  solve  the  equation  for  <p  numerically 
(or  analytically)  then  the  solution  will  obviously  become  very 
large.  A  noniterative  application  of  the  Klunker -Traci  procedure 
in  this  form  is  made  difficult  by  the  fact  that  one  must  first 
compute  the  flow  field  with  conditions  at  the  outer  boundary 
containing  the  unknown  surface  potential,  before  one  can  apply  the 
condition  that  the  potential  obtained  from  the  flow  field  computa¬ 
tion  must  equal  the  potential  assumed  in  assigning  the  far 
field  conditions.  For  noniterative  applications,  the  far  field 
condition  2  may  be  more  suitable. 

The  nonconvergence  of  the  iterative  process  can  probably 
be  avoided  if  one  derives  from  the  values  of  <j>  and  <|>  determined 
by  Green's  formula  (here  Eqs.  (70)  and  (71))  a  different 
boundary  condition.  For  the  present  problem,  this  possibility 
has  been  discussed  by  the  author  with  the  -iu<f>  +  <)>n  number  prescribed 
according  to  the  values  computed  from  Eqs.  (70)  and  (71) .  In 
this  formulation,  the  occurrence  of  the  imaginary  unit  is  important. 

One  finds  that  the  denominator  of  the  last  expression  will  not 
vanish.  This  would  probably  allow  one  to  use  an  iterative  procedure 
in  conjunction  with  the  Klunker-Traci  method.  However,  the 
difficulty  of  nonconvergence  still  exists  if  one  tries  to  carry 
out  the  flow  field  calculations  for  these  boundary  conditions  by 
an  iterative  process. 

For  a  case  with  spherical  boundaries,  the  situation  can  be 
discussed,  even  in  a  more  general  case,  by  means  of  the  particular 

solutions  (Eqs.  (33))  where  the  function  is  determined  by  Eq.  (28). 

2  2 

This  equation  shows  that  for  r  y  <  n(n  +  1)  in  the  three  dimensional 
case)  the  functions  are  nonoscillatory .  A  similar  result  applies, 
of  course,  for  the  two  dimensional  case  (ry  <  m) .  Eigensolutions 
become  possible  for  values  of  n  where  the  solutions  are  in 

the  oscillatory  region.  In  essence,  this  amounts  to 

r  fit  ■*  • ) 


35 


For  n  =  0  one  obtains 

)  <  *A 

It  follows  that  the  difficulty  encountered  by  the  presence  of 
eigensolutions  depends  in  essence  on  the  size  of  the  computed  part 
of  the  flow  field  (here  u  r^  or  u (r1  -  rQ)  ) ,  not  on  the  size  of 
the  body  (here  rQ) . 

Accordingly,  it  is  desirable  to  keep  the  computed  part  of 
the  flow  field  as  small  as  possible.  The  limiting  factor  is  the 
function  h(r)  in  Eq .  (54)  which  must  be  zero  at  the  location 
where  the  farfield  condition  is  applied  (except  in  the  case  of 
Klunker's  formulation).  In  realistic  cases,  the  function  depends 
also  upon  <|>  and  this  determines  the  size  of  the  flow  field  to  be 
computed . 
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SECTION  VII 


REMARKS  ABOUT  THE  NUMERICAL  SOLUTION  OF  THE 
PARTIAL  DIFFERENTIAL  EQUATION 

For  y  not  small  it  is  necessary  to  solve  the  flow  field 
by  direct  methods.  Eq.  (9)  shows  that  for  fixed  circular 
frequency  u,  y  becomes  very  large  as  the  Mach  number  approaches 
1.  Here  the  idea  of  Achi  Brand  (Ref.  3  )  may  be  useful.  One 
computes  the  flow  field  in  a  mesh  which  is  just  fine  enough  to 
include  all  eigenfunctions  which  are  responsible  for  the  failure 
of  an  iterative  process  and  uses  a  finer  mesh  to  derive  more 
accurate  solutions  by  iterations.  The  iterations  will  not  converge 
because  the  contributions  of  the  low  eigenvalues  are  not  completely 
removed  and  new  contributions  are  brought  in  by  the  iteration 
process  itself.  For  this  reason,  one  alternates  between  the 
determination  of  a  correction  by  means  of  a  direct  method  in  a  coarse 
mesh  and  an  improvement  of  the  solutions  by  iterations  in  a  fine 
mesh. 

Some  thought  should  be  given  to  the  choice  of  the  mesh. 

The  mesh  size  determines  the  waves  in  the  solutions  which  are 
picked  up  by  the  procedure.  The  mesh  size  is,  therefore, 
determined  by  the  value  of  y.  Without  further  information,  this 
gives  a  lower  limit  for  the  mesh  size.  This  limitation  applies 
to  the  entire  flow  field.  It  can  not  be  disregarded  in  distant 
parts  of  the  flow  field.  The  situation  is  different  in  steady 
problems . 

The  form  of  the  particular  solutions  (Eq.  (33))  shows 
that  the  main  variation  occurs  in  the  r  direction.  The  function 
Pmn(5)  cos (m6)  remains  the  same  along  rays  through  the  origin. 

This  means  that  these  particular  solutions  can  be  properly  represented 
in  a  grid  system  formed  by  the  intersection  of  equidistant  surfaces 
r  =  const  with  rays  through  the  origin. 

The  same  conclusion  is  obtained  by  the  following  consideration. 
Along  the  surface  the  potential  has  only  a  finite  phase  difference 
mainly  determined  by  the  boundary  conditions.  One  can  then  draw 
lines  of  approximately  equal  phase.  They  are  equidistant  curves 


to  the  surface.  For  a  plate  such  lines  are  shown  in  Figure  3. 

Along  each  of  these  lines  the  phase  differences  are  the  same. 

It  is  therefore  not  necessary  that  the  number  of  grid  points 
available  along  these  lines  increases  as  the  distance  from  the 
plate  increases.  The  number  of  grid  points  corresponds  to  the 
number  of  rays  through  the  origin  which  one  would  use  to 
approximate  the  particular  solutions  Eq.  (33)  or  their  two- 
dimensional  counterpart.  However,  the  distance  from  one  such 
line  to  the  next  is  critical  and  cannot  be  decreased  as  the  dis¬ 
tance  from  the  origin  increases.  Such  a  system  of  coordinates 
would,  for  instance,  be  approximated  by  a  system  of  confocal 
ellipses  and  confocal  hyperbola  as  shown  in  Figure  4.  Actually, 
the  ellipses  are  not  equidistant.  They  must  be  chosen  in  such 
a  manner  that  the  largest  distance  between  them  is  smaller  than 
the  critical  wave  length.  Writing  down  the  difference  equations 
for  such  a  system,  one  obtains  a  large  block  tridiagonal  matrix. 

Now  one  eliminates  the  unknown  values  of  the  potential  starting 
at  the  surface  and  proceeding  from  one  ellipse  to  the  next.  At 
the  end,  one  is  left  with  relations  between  the  unknown  values  of 
<j>  for  the  last  two  ellipses. 

In  the  elimination  process  one  may  have  some  pivoting  problems. 
They  would  arise,  for  instance,  if  for  one  of  the  ellipses,  the 
value  of  y  under  consideration  is  an  eigenvalue  of  the  problem 
with  the  normal  derivative  prescribed  at  the  body  and  zero  value 
of  the  function  prescribed  at  the  particular  ellipse.  The  ultimate 
aim  is  to  find  an  expression  which  represents  <{>  and  <t>  at  the 
outer  boundary  (or  equivalently  <j>  along  the  last  two  ellipses)  . 

This  is  the  information  which  one  obtains  from  integrating  the 
flow  equation  with  the  prescribed  boundary  conditions  at  the 
body.  This  information  is  combined  with  the  far  field  conditions 
most  conveniently  in  the  form  2.  This  is  a  global  condition  which 
connects  all  values  of  <J>  and  <J>n  at  the  outer  boundary.  Corresponding 
to  the  finite  mesh  size,  one  will  choose  only  a  limited  number  of 
values  n  and  m.  From  these  conditions  one  then  computes  the  values 
of  <j>  and  <J>n  along  the  last  two  ellipses  and  then  by  back-substitution 
for  the  whole  flow  field. 
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Figure  4.  Equi-Focal  Ellipses  and  Hyperbola.  (They  give  a 
regular  system  of  grid  point  with  some  similarity 
to  that  of  Figure  3.) 


It  must  be  conceded  that  the  mesh  advocated  here  is  not 
applicable  in  all  cases.  For  instance ,  in  the  transonic  flow 
problem  one  may  prefer  a  mesh  which  takes  its  orientation  from 
the  free  stream  direction.  Moreover,  in  the  transonic  flow 
problem  one  may  have  a  fairly  rapid  change  of  the  coefficients  of 
the  linearized  differential  equation  for  unsteady  flows  and  the 
mesh  must  be  fine  enough  to  reflect  these  variations.  These 
questions  should  be  explored  by  numerical  experiments. 
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SECTION  VIII 
FINE  GRID  ITERATIONS 

Assume  that  the  problem  allows  the  application  of  Brand's 
procedure.  Then  one  will  use  iterations  carried  out  in  a  fine 
mesh  to  improve  approximation  obtained  in  a  coarse  mesh.  It  is 
fortunate  that  for  the  fine  grid,  the  far  field  conditions  have 
approximately  a  local  character  although  they  are,  in  principle, 
global  conditions.  (This  means  that  they  give  relations  between 
all  the  values  of  4>  and  <J>n  at  the  outer  boundary.) 

We  use,  for  this  purpose,  the  third  formulation  of  the  far 
field  conditions,  letting  the  surface  S  coincide  with  the  outer 
boundary  of  the  flow  field.  It  is  assumed  that  the  computation 
with  a  coarse  mesh  has  suppressed  long  wave  perturbations  which 
are  detrimental  to  the  convergence  for  iterations  in  a  fine  mesh. 

In  the  third  formulation,  the  far  field  conditions  are 

expressed  by  means  of  a  density  function  f  (r')  ,  which  so  far  is 

unknown.  Assume  that  in  the  k  iteration  step  we  have  some 
(k) 

approximation  f  to  this  density  function.  Then  one  makes  the 
following  computations: 

1.  Compute  the  values  of  <J>  at  the  outer  pertaining  boundary 

for  this  function  by  means  of  Eq.  (59). 

2.  Compute  the  values  of  <|>n  at  the  outer  boundary  by 

means  of  Eq.  (59) .  In  this  computation  a  limiting  process, 

familiar  from  the  treatment  of  solutions  of  the  Laplace  equation 

by  means  of  an  integral  equation  method  is  required.  This 

function  is  denoted  as  d>  . 

n , outer 

3.  Compute  the  inner  flow  field  (iteratively)  by  means 
of  some  numerical  method  using  the  values  of  <f>  obtained  in  step 
1  as  boundary  values. 

4.  Determine  from  the  solution  obtained  in  step  three, 

the  values  of  d>  at  the  outer  boundary;  this  constitutes  d> 

n  J  n, inner 
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5.  Compare  the  result  of  step  2  with  that  of  step  4.  If 
the  difference  is  sufficiently  small,  then  the  far  field  conditions 
are  satisfied.  If  this  is  not  the  case,  then  one  carries  out 
step  6 . 


6 .  Set  f 


(k+1) 


=  f 


00 


^nouter  ~  <*>ninner^ 

This  correction  is  suggested  by  the  analogy  of  the  Laplace  equation. 
This  correction  to  the  pole  density  gives  exactly  the  jump  of  the 
normal  derivatives. 


The  procedure  is  effective  only  for  wave  lengths  which  are 
short  enough  so  that  the  corresponding  particular  solutions  (Eq. 

(33))  behave  in  essence  like  those  of  the  Laplace  equation.  We 
illustrate  the  procedure  for  the  Laplace  equation  in  two  dimensions. 
Let  us  assume  that  we  have  an  error  in  the  density  function  given  by 
6f^  =  c  •  sin  (m0) 

In  the  manner  discussed  in  Section  V,  one  obtains,  as  a  result  of 
the  integration  for  the  far  field  conditions,  a  function  6<J>  which 
is  free  of  singularities  inside  and  outside  of  the  circle  with 
radius  r^,  is  continuous  along  the  circle,  and  has  a  jump  of  the 
normal  derivative  given  by  the  function  6f.  Such  an  expression 
is  given  by 

ir  (r/r,)** t)  ;  r  <  r, 
r,  ( r,/r)  fis*{ ►*»  & )  j  r>y/ 

One  therefore  obtains,  as  a  result  of  step  1,  an  error  in  <J> 


n,  outer 


rt  &*(**&) 

and  as  a  result  of  step  2 ,  an  error  in  <t>r 

In  step  3,  one  solves  the  inner  problem  with  the  conditions 
6<f>  =0  for  r  =  rQ,  and  6 <f>  is  given  by  the  result  of  step  1: 

(%r+ 


r  j.  =  —  'll/*  fa*  &) 


rf  /**/**  &} 
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This  then  gives  for  the  normal  derivative  at  r  =  r1 
***** 

and  for  the  difference  between  the  normal  derivatives  of  the  outer 
and  the  inner  solution 


,  >-(%i  . 

*>[' '  77757-/'  <■*»(”*>  ,tfx 

I  f!  '  U  * 


y 


j»i 


The  correction  carried  out  in  step  6,  then  leads  to  an  error  in  the 
density  formula  given  by 

y  ^ 


This  amounts,  indeed,  to  a  considerable  reduction  of  the  error. 


SECTION  IX 
CONCLUSIONS 

The  results  needed  for  the  practical  application  are  given 
in  Section  V.  The  Klunker- Tracy  iteration  is  useful  for  small 
frequencies,  but  other  methods,  for  instance  method  3,  may  be 
equally  efficient.  Most  useful  is  probably  a  combination  of  the 
formulations  2  and  3.  The  applications  are  demonstrated  in  an 
overly  simple  example  in  Section  VI.  This  example  may  not  be  very 
useful  for  practical  work,  but  it  shows  possible  difficulties,  and 
properly  interpreted,  it  also  shows  the  steps  that  are  to  be  taken. 
Section  VII  tries  to  counteract  the  strong  analytical  flavor  of 
Section  VI  by  discussing  in  general  terms  which  steps  are  to  be 
taken  in  the  frame  work  of  a  finite  difference  method.  The 
author  believes  that  the  ideas  of  Achi  Brand  who  advocates  the 
simultaneous  use  of  grid  systems  of  different  size  can  be  very 
useful  in  practical  work.  Section  VIII  describes  the  use  of 
the  third  formulation  of  the  far  field  conditions  in  an  iterative 
solution  for  the  fine  grid. 

These  discussions  are  carried  out  with  a  view  toward 
practical  numerical  application.  The  step  to  a  workable  program 
in  a  specific  situation  has  not  been  carried  out  so  far. 
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APPENDIX 

RELATIONS  CONNECTING  DIFFERENT  PARTICULAR  SOLUTIONS 


The  functions  f^  are  defined  by  Eqs. 
One  has 

,(H>)  -ft-/  .  r-1  -H.  A 

/  (r)  c  T  Jtx/3  (-  i/±r )  2.,  O',  r 

6  k*o  * 

with 


/  yv!  ( lK-k )! 

k,  ■'  (n-k)!  fen,)! 


(29),  (30) 


Notice  that  for  all  n 


and  (32). 

(Al) 

(A2 ) 


\  J 


a.^  s  2.^/u. 


The  inllowing  relations  will  be  derived 
d /  ,  JntO  A  **" 

-  -(*'  +  /){  +  7~TT  / 


(fr 


4-h*--/ 


4rr-  / 


M+O 


jUL 


/ 


(A3) 


(A4) 

(A5) 


This  is  done  by  substituting  the  power  series  development  into 
these  formulae.  Notice  that  one  has  two  contributions  to  df^/dr 
one  from  differentiating  the  power  series,  the  other  from 
differentiating  the  exponential  function.  One  obtains  for  the 
coefficient  of  the  term  with  the  power  r-11'1  +  k  of  dfn/dr 


k! 

fn,-k  -O- 


(lK-k-/)! 


(**)! 


f-Ktk,) 


/  */! 

£7  fn-k)f 


fen,  -k ).' ) 
fen,)!  / 


-  fefe'f"  '  -  f-n?-  *'+'*'J*  ~  ~T 

1  f  (*,+  >)!  fo-V!  (in,)'  ( 


The  coefficient  of  the  term  with  the  power  r~n  of  dfn/dr  + 
n4*  1 

(n+l)  f  is  computed  next.  One  obtains 


UM)!  fo-k)! 


U-H)!  fr-k)!  (L*,+l) 


The  term  within  the  braces  is 


The  term  within  the  bracket  vanishes  for  k  =  -1  and  k  =  0  because 
of  Eq.  (A3) .  It  is  a  quadratic  function  in  k.  Therefore,  it 
can  be  written  as  k(k+l)  time  constant.  The  constant  may  depend 
upon  n.  It  is  readily  obtained  by  letting  k  tend  to  infinity. 

One  obtains 

f  t  x  ~  k  (6,+t ) 


The  coefficient  of  the  term  with  power  r 
fn+^  is  therefore  given  by 

The  coefficient  of  the  terms  with  power 
n-1 
ak-l 


of  df  /dr  +  (n+1) 


comes  from 


This  suffices  to  establish  Eq.  (A4).  Eq.  (A5)  is  verified  in  an 
analogous  manner.  The  coefficient  of  the  term  with  power 
r  n  2+k  of  fn+1  -  (i/r)fn  is  given  by 

*  a  * ‘  t fa-i+t)!  (ik+i)!  fo-AJ/  fa*,).1  * 

.  At  '*>/  (i*, -l)/  , 


■> 

The  term  in  braces  vanishes  for  k  =  0  and  k  =  1  and  behaves  as  k“ 
for  k  -*■  00 .  Hence, 


and 


A- 


H,-#/ 

A. 


/Jm  ) 

4  6  (i_d!  {„ -in)!  fe+o+n!  * 


The  coefficient  with  the  power  r  n  of  fn  ^  is  given  by 

n-1 
ak-2 


h/~> 

Co 

1-2 


!  /  (H>  ->)!  fen-i.)! 

fa-k,*/)!  fen-2-)! 


From  the  last  two  formulae,  Eq.  (A5)  is  immediately  obtained. 

From  a  theoretical  point  of  view  the  existence  of  formulae  of  the 
type  A4  and  A5  is  not  surprising.  One  deals  with  contiguous 
confluent  hypergeometric  functions. 


Similar  relations  exist  for  the  Legendre  functions.  We  use 
the  definitions  given  in  Reference  4 ,  but  take  the  formulae  from 
Reference  5.  There  exist  a  number  of  different  definitions  in  the 
literature  (also,  there  is  probably  a  misprint  in  the  section  called 
notations  in  Chapter  8  of  Reference  5).  For  clarity,  the  definitions 
to  be  used  and  the  basic  formulae  are  repeated  here. 

ft).Z  _  (A6) 

r*0  l^r!  (»-r)!  fn-2--r)! 

Reference  4,  page  274.  We  have  not  given  an  upper  limit  for  r. 

The  sum  terminates  automatically  when  one  of  the  factorials  in 
the  denominator  assumes  a  negative  argument  (here  for  r  =  ^  +  1 
if  n  is  even;  and  for  r  =  — if  n  is  odd)  .  That  is,  the  sum  goes 
up  to  r  =  j  for  even  n,  and  up  to  r  =  for  odd  n. 

Furthermore , 


I 


tto  '“•'A  d  ?  (b) 

l  *  (-  f*  ft-  ?)  - tudLL 


Differences  in  definition  may  occur  because  one  may  replace  (-1) m 
by  1,  and  because  of  the  argument  (l-£2) (in  the  formulae  of 
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m 


Reference  4  one  has  always  (c  -1)).  This  would  introduce  a  factor 
(+i)m.  For  this  reason,  the  formulae  of  Reference  5,  made  con¬ 
sistent  with  the  present  definitions,  are  repeated  here.  The 
signs  have  been  checked  by  evaluating  these  expressions  for  the 
highest  powers  of  £.  Also,  the  lowest  powers  have  been 
examined . 


2^'(i )  -  £a)  -  r*+**)£ s&> 


■'A 


(A8 ) 


-—y* - /  -  (k  +  *»){*-»,+>)  fr'fj 

*  S 


from  (8.5.1)  of  Reference  5 


(A9) 


from  (8.5.2)  of  Reference  5 


-L~  Uk-k-h)?*' 


(A10) 


from  (8.5.3)  of  Reference  5 

(no  check  necessary  since  the  superscript  of  Pm  does  not  change) 


(f-f) 


m  42+ 


A 


-  2  *  (*, +  *»)?** 


(All) 

from  (8.5.4)  of  Reference  5 
(no  check  necessary,  since  the  superscript  of  P™  does  not  change) 


r  -f  a 

Hrf/  •*-! 


(A12 ) 


from  (8.5.5)  of  Reference  5 
The  following  formulae  will  be  used. 


(A13) 


<!£ 


*fr*or***’>£] 


(A14) 


(•-&*+  -  a 


■*  /  **+/  m-t  •* 


(A15) 
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j?-  -zrnf-M' S') 


/*// 


0~£)  -J {'*>-*+ 1 y*-n,+2-)?4i't.  -  fa +»v -/]{*,+»*) t  J 


(A16) 


(A17) 


S  ^77  ffa'*+')fo-*+*)[+f''fJ++'7  £7 

+(*v  tm  -t)fi  /Ht )  J 


(A18) 


Eq.  (A13)  is  identical  with  Eq.  (A10) .  Eq.  (A14)  follows  from 
Eq.  (All)  with  some  further  computations  by  eliminating  c;P™(s) 
by  means  of  Eq.  (A10) .  Eq.  (A15)  is  identical  with  Eq.  (A12) , 
with  m  replaced  by  m  +  1 .  Eq.  (A16)  requires  some  further 
arimetic.  One  obtains  from  Eq.  (A7) 

d  i  =  Trp-S*  ~  (/-fj'k  *y 

Hence , 

d  S 

Now  from  Eq.  (A13)  with  m  replaced  by  m  +  1 

-  ; - /{**■-**')?  +  J 

°  **  Z»-n  L  *  +t  m  ' 

Furthermore,  we  use  Eq.  (A15) .  Substituting  these  two  expressions 
into  the  last  equation,  we  ultimately  obtain  Eq.  (A16) .  To  obtain 
Eq.  (A17)  one  first  writes  down  Eq.  (A8)  with  m  replaced  by  m-1 

('-iVAC  ’!(*■—*•> STZ"(y  - 

and  substitutes  Eq.  (A10)  with  m  replaced  by  m-1 


«•"  Z»+/ 


- — //**-** -*2)  •?  V  (*t  +  *n-ij T  J 
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To  derive  Eq.  (A18)  we  start  from  Eq.  (A9) 

(/-£)  *  (»+**)(* -■*•+ 00 -&*)'*’* 
Now  from  Eq.  (A8)  with  m  replaced  by  m-1 


Now  from  Eq.  (A10)  with  m  replaced  by  m-1 
*  r'"‘  * (***->£. 


This  leads  to  Eq.  (A18) . 

These  equations  are  now  used  to  derive  formulae  which  express 
x,  y,  or  z  derivatives  of  certain  particular  solutions  by  other 
particular  solutions  with  changed  m  or  n.  We  have 


7>i  m  j_  _ 
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■  p  .v  ja«* 


(A20) 

(concluded) 


Using  the  formulae  just  derived,  one  obtains 


yj-  ^ 


Here  Eqs.  (A4) ,  (A5) ,  (A13)  and  (A14)  are  substituted. 

jc*i  +, y  I  fr  '»('»-», +!) 1  f *>+!){»+**) y*  If  f  /*~7 

l I [}  *»■>>  *-*Ji  +n*-j  r 


t 


ft*-  **+  t  w  ■>> y->n 


Here  the  products  P™+1  fn_1  and  P™_1  fn+1  vanish  and  one  finds 


fe  vyv/ZZV^  */-  ",  - 
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One  observes  that  =  0  if  n  <  m;  that  is,  this  sequence  of 
equations  starts  with  m  =  n  and  then  expresses  higher  derivatives 
in  terms  of  function  P™  with  higher  values  of  n. 

Next  consider 

«* 

4  CtttHfjM* 


Hence , 


£  K(Vf  -  37/ "** 
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(A22 ) 


2-[T~(!Jf''''r>c«(~+il+ 

-  *  ~  tZdto-^J^ 

t  tZivc-i1^  £’}  ^">t> 


(A23) 


Now  we  substitute  into  Eq.  (A22)  the  expressions  Eqs.  (A15)  and 
(A16)  and  (A4)  and  (A5) . 
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Hence 


4- 


(A24) 
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Take  the  expression  (A23)  and  substitute  Eqs.  (A17)  and  (A18) 

j;  [K(Vf™r)  ("f*] 

y/^y-*!  -/ )(*v  +  **< )  [ fa-t/J  "  +>*/ 

y  jzji  l(* - *  -*• * ^ 

+  [(*-»*')(*  +  -  fat  m -!)(•*+»+)?* f]  x. 


*■  ■*•»+/  /'fcvy 
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«  -7— — ^ jf-  !)(*,- 2  -  fa *»){*,+ /)?  /  * 

*»»■/  /  tL  *¥/  '  <*-/'' 

//  (r)l 

4 [ (*-**  +  !){'*'- **  -(*  +  "*-/)(*  +  **')  _/  /  * 

r  /  ,  rf*’t,l  .  ft1*-  ,  7) 

[-(*  +  !)  f  rr>+^,f  MJJ 

tn-l  fall)  %-/  .^w/ 

»  eoi((»-t)y)  [- +  ,  /  -  /«,+m-t)fa+M,)2"  ^  / 


Hence 


-^r  llZWf'^7)  coifay-)]*  jj  [ )] 

•H-/  f**>  ,  y*-^  7 

e  +  i)(*i'i»  +L)I  +/t)f  -fa*1*  -t  )/''»•***’)  r  f*Jj  * 

X  Coi((**-l)p  )  (A25) 


There  exist,  of  course,  analogous  formulae  with  sine  and  cosine 
interchanged  (and  certain  changes  of  signs) . 

Similar,  but  much  simpler  relations  exist  for  the  two- 
dimensional  case.  One  has 


x  *•-  xA*^ 


,  <9  »  axct${y/x) 

Gos  ft  j  -7-  j*  St*t  ft 

77  -  -  -pV^  ,  yCoj  ft 


7T 

Vft 


Then  one  finds 
?  r  ><  M. 


n  r  , ,  f *■/  7  ^ 

7*  ^  f^)  cot  fa  ft)]  ^  ^ifafticosft 


■*  r  ^  r)jcy>/»t  ftju  n  ft 
Cflf**  ft )  tol  ft 
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Hence 


*)] "  7y[ r)  /*  *)] 

,  ,.fo  .  . 


f*T 


(A26 ) 


J7 1  r)  Ctxfa  6)]  +  3j  [#£(#')  *»  (*"  &j] 

,  fu,  // 'V T)]<u>i[fa-,J*]  (A27) 

Lc  T  **'  “ 

Consider  some  linear  combination 

2  fax)  -  5  «£  (•f'X)  +  c^/VpUx) 

(c1  and  c2  constant  real  or  complex  coefficients,  J  and  N  are 
respectively  Bessel  and  Neumann  functions  of  order  p) .  One  has 
the  following  relations 

■£ /*  *•$&*»  *  *  **$.,  '<■*>  -f  $M»)  *  * 

•-**'$,/**>  -%$/**>  +  * 

r 


fox) 


a'x 


Hm  is  such  a  linear  combination  of  Jm  and  n  .  One  therefore 
obtains  from  Eqs.  (A26)  and  (A27) 


Ji  iH^t^rlcahulf)]-  jL **.(*, !■)] *  un^MtOt-]  (A28) 

■k iCr)u,fm^  (a29| 


There  exist,  of  course,  analogous  formulae  with  sine  and  cosine 
interchanged  (and  some  changes  in  sign) . 
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